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1. INTRODUCTION 

Electric Power Research Institute (EPRI) carried out a large survey on faults in induction motor 
in1985. Bearing faults contain over 40% of all machine failure. For this reason, the detection of bearing faults 
has been an important research area [1-6].The different faults occurring in a rolling-element bearing can be 
classified according to the damaged element as [7]: 


a. Outer raceway defect 
b. Inner raceway defect 
c. Ball defect 

d. Cage defect 


Bearing faults as all other induction motor failures ought to be diagnosing as early as possible in 
order to prevent the complete failure of the system and unexpected production costs. For that, the challenge is 
to have an adequate model of the induction motor which behaves as the real motor, and is able to generate the 
desired signals as we see in real motors under healthy or faulty state [8], [9], or in motor-pump system [10]. 
The model for bearing fault proposed by Schoen has been applied in several works as, [11], [12]. Blödt [13] 
consider this model incomplete, and extended it by the consideration of the torque variations, and the 
introduction of the radial movement of the rotor center. The fault impact on the airgap length is considered 
by a series of Dirac generalized functions of vibrations analysis and has been used in Blédt’s model [13], 
[14]. In this approach, one can notice that, firstly, the surface of the Dirac function is equal to 1, but its 
amplitude tends towards the infinite one, and that makes the simulation of inductances impossible, and 
secondly, it’s true that the contact angle is considered, but the angle of the defect does not exist. Then, for 
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these two raisons, we consider that this approach is not highly reliable. In this work a new approach with 
taking into account the angle of the defect will be proposed. 

In this paper, one of the fourth types of bearing faults classified above is studied, namely the inner 
race bearing fault, and because the similarity which exists between this last and the dynamic eccentricity 
fault, a new way for formulating a general model is suggested. The main objective is to offer accurate 
mathematical formulas of all the inductances of the machine in order to simulate efficiently the behavior of 
the motor under either healthy and/or faulty conditions. 

Bearing faults can be diagnosed with help of many condition monitoring methods such as vibration 
monitoring [15], [16]. This method is applied with the use of broad band, narrow band and the spectral 
analysis of the vibration energy of motor. But this method has a drawback of high cost as it requires 
expensive accelerometers.Thermal monitoring [17]. This method is done by measuring the local and the bulk 
temperature of motors; it also can be achieved by parameter estimation. This method has a drawback that it 
might be too slow to detect fault. Current monitoring [18-21]. The method used for diagnosis [22] is accurate 
and remains one of the most used, but require expensive sensors or specialized tools. For that, and because 
the stator current harmonics measurement can provide the same indications without requiring access to the 
motor, in this paper, stator current monitoring is used. 


2. DYNAMIC ECCENTRICITY AND INNER RACE BEARING FAULT 

In (DE), the center of the rotor is not at the centre of the rotation and the position of the minimum 
air-gap rotates with the rotor (Fig1). 

With the same manner used in. The 3D expression of the air gap variation is [23]: 


Bionic 0, ? x) = Eo 1 Ô; 1 ` cos(@ 6,) a) 
i lne) 


And from the approximation for the first harmonic, one can get: 


g3. (0,0,,x)= A i +6, i = =) cos(-6, ) (2) 


Eo (DE) 


Figure 1 shows the two types of dynamic eccentricity (uniform and non-uniform) 


Figure 1. Non uniform (left) and uniform (right) dynamic air gap eccentricity 


dq: is the relative degree of eccentricity. 
The variation of 5g versus the axial position x is: 


PEP 3) 


The equivalent level of dynamic eccentricity is defined as a function of the axial position x by: 


5,(0,0,,x)= afi- 7 Jono 6.) (4) 


ex 


Analytical Model of Cage Induction Machine Dedicated to the Study of the Inner Race ... (Hadjami M) 


460 o ISSN: 2088-8708 
If the variation of “dq” versus the axial position x is removed, we return to a 2D expression. 


2.1. Inner Race bearing Fault 

Bearing is an arrangement of two concentric rings with a set of balls which spin in raceways 
between the inner ring and outer ring as shown in Figure 2. When a bearing spins, any irregularity in the 
inner raceway surface generates vibration occuring at periodic frequency, called (BPFI) Ball Pass Frequency 
of the inner Race, and is created when all the rolling elements roll across a defect in the inner race as shown 
in Figure 3. 


Beg p vus Outer race 


K bldh 


Inner race 


Figure 2. Bearing geometery Figure 3. Dynamic eccentricity caused by the inner 
raceway defect 


6,(0,0,.x, fy) is the equivalent temporary level of eccentricity and y is the defect angle. 


Cages function to maintain all the balls at a uniform pitch and rotates at the fundamental cage frequency 
cage, Which is given by [24] 


1 D 
Orage =~ @,| 1+— cos (5) 
a51 Beco] 
B: is the contact angle. 

Defect frequency depends upon the bearing geometry and shaft speed @,. In particular for the inner race, and 
running at a certain rotor speed œw,, the predictable characteristic fault frequency fint is [24]: 


fa = sis Pecos) 1) 


Npa: The number of balls. 

The contact between the balls and the defect caracterised by the angle y, generates a temporary static 
eccentricity. The new way for modelling this fault impact, is the introduction of a digital square wave 
Figure 4. This signal is expressed as complex Fourier series development. 


D„= an : >. : inl “Joo on 1) (7) 
~ an aaa 2 2 


Figure 4. Digital square wave signal 


Equations (4) to Equation (7) can be combined into a simplifie expression of the equivalent level of 
the outer race bearing fault. 
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8 (0,0 x)= 4 6,404 1-Ž S sin{ AZ Jos 0 0+4 Oy F z) (8) 
27 m l a 2 2 


ex J A=1 


By analogy with (1) and (2), the new formula of the air gap lenght, and the inverse of the air gap 
function approximated for a small variations, become respectively 


/ 6 al. 
8 (Dyn-ece (9,9, 5 x)= «I = 6, u: - Edsn( 22 eosa 6.4 (2=2)] (9) 


= 1 yo 6 x1. y y (10) 
! 0,0,,x)=—|1+-— 8, +—+|1 A 0-0,+4 0. 
eieaa Yt a Ahi- Ei 22 noo aa 


y = 0 (Healthy state), y = 2r (Dynamic eccentricity), y # 0,27 (Inner race bearing fault) 


3. INDUCTANCES CALCULATION 
According to winding function theory, the mutual inductance between any two windings i and j in 
any electric machine can be developed by [25-27]: 


2al 


„(0,0,,x)= Mor Tae ‚(0.0,,x)g" (0,0,, x)dxd0 (11) 


2a1 


f f n,(0,0,.x)¢ '(0,0,,.x)dxd0 


0,0,,x)=n,(6,0,,x)—22 (12) 
M*(0,6,,x)=n,(0,6,,x) rile 0-0) 
1 2al 
(¢'(0,6,,x) Paate (0,0,,x)dxd 0 (13) 


where 6, is the angular position of the rotor with respect to the stator frame, @ is a particular angular 
position along the stator inner surface,g“(0,0,,x) is termed the inverse air gap function. The terms n(0,0r,x) 
M(0,0r,x) are the turn function and the modified winding function in windings i and j respectively.(the 


modified winding function theory is used because the non uniformity of the air gap). The model equations of 
the induction machine are the same as [28]. 


3.1. Turn and Winding Rotor Functions 


3.1.1. Healthy Conditions 


In this case the air gap length g is uniform and equal to gp (6,=0), the turn function, and the winding 
function of the rotor loop k (with taking into acount the skew angle 6, of the rotor bars) are respectively: 


1, (0,0,,X) = alae z sinh Z hos h| 6-86. [x e+% x (14) 
22 Aan 2 2 l 


a 2&1./( a 1 
n (0,0, x)= F i 25 isna ISLC 0, (x A n <) (15) 


a,: The rotor loop opening. 


2m 


n, 
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np: The number of rotor bars. 
The variation in the angle y, versus the axial position x i 


ô (paces (17) 


3.1.2. Inner Race bearing Fault 
The modified rotor winding functionM,.;(0, 6,, x) can be expressed independently from the value of 


P as: 


sin| 2 
SE a $) col) l ao {e-5)e 4) n 


3.2. Turn and Winding Stator Functions 
3.2.1. Healthy Conditions 
As the same as for the rotor (6d=0), the turn and the winding functions are respectively [28]: 


Na, 2N,ck 27 
n (9)= ++ += “ht scos| hp| 0-0, -lq -1)— (19) 
O) Ng MSE scr (0-04 —(a-) | 
Na, 2N,ck,, 27 20) 
0 1i } t wh h 0-0 = =f panies ( 
Neal ) on EA h se A| 0 (q =) 
21 
fa airoja a, “Fe. N, = NN, Ds Kun =Kankins Kan =sin( 22 J) vn( ME in "2E } a,= z } e 


were, a,: the stator coil opening, Q: the nuber of slots per pitch turn, N,: the number of the conductor per 
stator slot, Ne: the number of slots per pole and per phase, N;: tne number of stator turns in series, and p: is 
the number of pole pairs. 


3.2.2. Inner race bearing Fault 
One can notice that the modified winding function M,, (O) is different from the winding function 


N,q(9) only in the case of p=1. 


M sp l0)=N,(0) (22) 


2prl N?82k? iy 2a) fy 2&1. (7 7v\]| 23) 
Negin)(9) = N (0) > - O1 a, | [oO % l-07 Sigs 7 sinl AS J008} Al An 5 


3.3. Derivation of Inductances 
3.3.1. Healthy Conditions 
Using relation (12) and when the skewing effect is considered, the different inductances have 


similar formulas as presented without skewing effect [28], except that the stator-rotor mutual inductances are 


multiplied by a skewing factor kins for each harmonic component h: 


Kilt, = sin hp 2) / (w 2) (24) 
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1 2 
ARET = 2M; os co 0,- -O,+ i(k - 5 AG) z) (25) 


with [ur z sag l S = el 
0 


3.3.2. Inner race bearing Fault: 
By using (14) one can obtain the different formulas of the stator and rotor inductances. The stator 
and the stator-rotor mutual inductances are dependent on the value of p, (p=1 or p1). 


2url 1 y = | [2 z) í z) 
M (pe) = sin cos} A} Ou -= 
rjk(DE) g n, a 2 in 2 


<{afemla, 0.) coslo, ~@)anle,~€))-2f esl, -0,)+ Zeosl,~@ ose, -0)}}29 


with: 

1 a 27 1 ki Ko 

O ra 0, H: =} 2 3 Pig A, A, (q 1) 3p k Osj a (A + 2 (j + k ler, 2 

k 
me ST i asain fei fs. B= [sil Jeo Z) -1 
i & 0 T hp 2l 2 JA r Le 2 2 
Lr jk(healthy)> Ls, 192(heatthy) ate the inductances given in the case of symmetrical rotor respectively 

[28]. 


a. The case of p#/: 


Lsq1a2 = Lsa1d2(healthy) (27) 


sr sr Ou lL S ayer ske 2 Pile s 1 
Lorann) = Estgi(neatiy e 7 fi- D porien (to a%-(4 D ) F (np +(-1) |: e 
ô, Y l — Mis 1 ) (ip +(-1¥}, ske “ Qa Y, 5 A 1 
K 2 z pare iy [oo 2 Kiti) x sin} hp a Bay = 5 = (hp +( 1) [«-3}e 
16, Lae. 2 = Fpi 
a a T E E O G ar |i-g 
16, l > Dsl 22) M p(y) cos (np +(- 17 y, kj 
m 2| Vd Jie A 2) hp+(-1) > (iy) 


sl pega a Ezoza kalpat e-i] ] (28) 


b. The case of p=/: 


For this case, one can use the inductances aR 1Sq2(Dyn-eco) ° Lkahpyn- 2 developed in our old works 


cited previously, with the same manner as for the case of p+1 by adding +0, . 
The stator and rotor self magnetizing inductances are respectively L,,,,; (=, 2, 3) and Lj G=1: n,-1). 


The different inductances formulas presented for the case of the Inner race bearing fault confirme that, when: 


y=0 => Lrjk = Ly jk(healthy) (Healthy state) [28] 
Y=20 => Ly je = Lrjk(Dyn-ecc) > Lkqh = Lkqahpyn-eco) (Dynamic eccentricity) [29] without skew 
y 0,27 => Lrjk = Ly jxcinty > Lkqgh = Leanne) (Inner race bearing fault) 


Figure 5 shows that when the rolling elements roll across a defect in the inner race, a dynamic 
eccentricity is created at periodic frequency called (BPFI). 
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Figure 6 shows the effect of the defect angle on the eccentricity created by the inner race bearing 
fault on the mutual inductances, and one can notice that the enlarging of the defect angle, involves an 
increase in the duration of the fault, and if the value of the defect angle is very small (close to zero), the 
variation in the inductance waveform is almost invisible. 


— Inner fault 


~~~ healthy 
~~~ Dyn-ecc 


Figure = Mrama a Lon au dynamic Figure 6. Mutual inductances L,),; with y and 6d =0.3 
variationseccentricity (red) and inner race bearing a=) 


fault (blue) 


Finally, the Figure 5 and Figure 6 show that inner race bearing fault causes an asymmetrical mutual 
inductance between stator phases and rotor loops. 


4. INNER RACE BEARING FAULT SIGNATURES IN STATOR CURRENT 

To find in the current of the stator the different signatures of the inner raceway bearing fault, and to 
understand their origins, it is first necessary to identify the harmonics induced by this defect in the rotor 
current.For an ideal three-phase supply, the stator’s currents which flow in the stator phases are of the form: 


[7 l= Lon of os — (q -1)) (29) 


The derivative of rotor flux in vector matrix form can be written as: 


4 l=“ TAIA EE TAA (30) 


d” a 


Where: 
@, : is the main voltage supply pulsation 
[L] : is the rotor loop inductances matrix 


[uz | : is the stator to rotor mutual inductances matrix 


Since the skewing of the rotor bars does not generate any new harmonics, so the study will be made 
with the mutual stator-rotor inductance without skewing effect, and which can be expressed as follows: 


P F y ô, L Aini 2a \_ : 1 
Lg = Baty + 2 EIS cys (a ina ¥ (hp +( 1) k 2 a, 


af E Cae 0, -(q 0%} 10,422 ¥ (hp +( mG sje) (31) 


m2 lex) h=l A=1 


After development of rotor flux derivative and some simplifications we get: 


[z] = rant Lean + aca) | (32) 


where: 
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le pls 2i cos(s,, ene ke, -hpo) (33) 


k=0,1,..n,-1 


[La y= vr cos(s ska, —hpg) (34) 


h=1 
k=0,1,..n,-1 


lal- X cos(syy,t + (hp -1)F2pka, ¥ A2af,,, —hpg) (35) 


h=1 
k=0,1,..n,—1 


All the harmonics induced in the rotor current have their images at the stator current. From there, 
and following the same previous steps, we can write the derivative of the stator flux as follows: 


lv.) L] A) + Mal, deal alal G6) 


were: 


[ku] ]_ alten lEn p] dlya] ]_ dlri nal (‘)) Hepa) Helle (37) 
dt dt dt dt dt dt 


After development, we find that the terms in (40) are always null except for the follow harmonic 
orders [29], [30]: 


[i= h -[ 2+) , h {2 i a] „h (2 = v1] | (38) 
P A=I,2,.. P A=l,2,.. P Aih 


The voltages induced in the stator produce a new series of high frequency in the stator current, and 
the expression of these components are: 


a (220-21) fF a f, Fkfa A=1,2,3,..k =1,2,3,. (39) 
p p 


The voltages induced in the stator produce a new series of high frequency in the stator current, and 
the expression of these components are: 


faa = (A s)+ ] fee) BAIS AEN (40) 
p p 


The Table 1 shows the new expression of the frequencies found by our approach, and that found by 
the approach of Blodt. 


Table 1. Stator Current Inner Race Bearing Fault Frequency 
Inner race bearing fault frequency 


The Blödt approach NEE = Í, + J + kf nt 


The present approach ames = a (1 — s)F if F 1 Í, F kf 
P 
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5. TEST LABORATORY 

A test bench have been set up at Ampere Laboratory CNRS,UMR 5005 of Claude Bernard 
university Lyon 1, France, in order to validate the previously discussed analytical model. The Tested machine 
is a star connected induction motor. Ratings of the tested motor are a 5.5kW three-phase induction machine 
with 48 stator slots, 28 rotor bars, and 2 pair poles from 5 Leroy Somer manufacturer,LS 132S with a 
protection factor of 55 (protection against dust and against water jets), nominal motor speed is 1440tr / min 
for 11.4A and cosp=0.84. Stator phase resistance 1.315Q. It has two NUP 206 type rolling ball bearings with 
13 balls. 

The defect was simulated on a roller bearing by the addition of a contaminant (cement powder). The 
inner ring fault (a stripe on the inner ring) was carried out on 3 levels of severity as shown in Figure 7. 


4 


Figure 7. Failure on the inner ring of the bearing NUP 206 (a-Level 1/3, b- Level 2/3 
and c- Level 3/3) 


For a load speed of 1406 tr / min for a level (1/3) of defective bearing, 1423 tr/min for the level (2/3) 
„and 1420 tr/min for the levels (3/3) and for 1420 tr/min the safe bearing - Bearing NUP 206, and using 
equation (53), it sets the table for the frequencies detected a small angle default (y) so a hole, or a larger 
defect angle therefore for our case a groove in the inner ring. 

Indeed, for any motor we always distinguish visible small eccentricity rotor frequencies, which 
explain the existence of some specific harmonics of the inner race bearing fault even for a healthy motor as 
shown in Figures 8 (a), (b), (c), Figures 9 (a), (b), (c), when the motor is tested with a defect angle at level 
(1/3). 


Aerie £ 
? 

$ 

i 


too tte ‘= 139 s. tee too tre tao oe tto tes mo ers ‘on ro irs +e ten 
Froauenoy (tes 


Figure 8. (a) Experimental FFT spectrum of the stator current under healthy condition (Blue), and faulty 
condition (Red), s=0:061. Level (1/3) between 100Hz and 180Hz (Zoom on (b) 113.03 Hz, (c) 172.21 Hz) 


Mota £ 


Seeeen ea: 


i 


Figure 9. (a) Experimental FFT spectrum of the stator current under healthy condition (Blue), and faulty 
condition (Red), s=0:061. Level (1/3) between 170 and 250 (Zoom on (b) 189.56 Hz, (c) 213.03 Hz) 
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This group of harmonics is summarized in Table 2. 


Table 2. Specific Harmonics of the Inner Race Bearing Fault in Stator Current (Level 1/3) 
AN 1 
Fintece -( Pe if Kha —> Fr 


P 
Healthy Faulty 
A=1,k=3 173.9 Hz 172.21 Hz 
114.43 Hz 113.03 Hz 
1=3, k=10 214.43 Hz 213.03 Hz 
190.76 Hz 189.56 Hz 


The spectrum of the faulty condition case (Level 1/3) contains new harmonic components (162.30 
Hz, 262.30 Hz, 202.83 Hz and 302.83 Hz) as illustrated in Figures 10 (a), (b), (c) and Figures 11 (a), (b), (c). 


SERJI 


Figure 10. (a) Experimental FFT spectrum of the stator current under healthy condition (Blue), and faulty 
condition (Red) s=0.063. Level (1/3) between 140Hz and 220Hz (Zoom on (b) 159.28 Hz, (c) 198.48 Hz) 


BARKA 


tezas 


s.. son see are san r.a tan 8 ton roe re 8 pvo fee poo sos 


Præuenoy te) 


Figure 11. (a) Experimental FFT spectrum of the stator current under healthy condition (Blue), and faulty 
condition (Red) s=0.063. Level (1/3) between 230Hz and 310Hz (Zoom on (b) 259.28 Hz, (c) 298.48 Hz) 


This group of harmonics is summarized in Table 3. 


Table 3. New Specific Harmonics of the Inner Race Bearing Fault in Stator Current (Level 1/3) 
AN 1 
inzece -( = (l-s) ip Fin È f. 

Fint-ece = 8A -8)-D f, — 6 fin, 1/2 f, =159.28Hz 
faee = 280-5) +) f, -6f,,, +1/2 f, = 259.28Hz 
Fimt-cce = 5601-5) -1) f, -13 fa -1/2 f, =198.48 Hz 
faee = 561-5) +1) f, -13f,,, -1/2 f. =298.48Az 


A=2, k=6 


A=4, k=13 
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The motor is now tested with a few significant increase in the degree of severity (level 2/3) from the 
previous (level 1/3), but with a clear reduction of the load relative to the previous test (1423 Tr/min instead of 
1406 Tr/min). 

Note that the amplitudes of the specific harmonics of the fault inner ring have decreased due to the 
reduction of the load, against the minimal increase in severity had no impact on either the amplitude of the 
harmonic and on the new generation of harmonics. The Figures 12 (a), (b), (c) show an example of a pair of 
harmonics (159.28 Hz and 259.28 Hz). 


testa a 
è 


Figure 12. (a) Experimental FFT spectrum of the stator current under Healthy condition (Black), faulty 
condition Level (1/3) (Blue)s =0.063, and Level (2/3) (red)s=0.051 between 150Hz and 280Hz (Zoom on 
(b) 159.28 Hz- 259.28 Hz Level (1/3) and (c) 162.73 Hz-262.73 Hz Level (2/3)) 


This pair of harmonics is summarized in Table 4. 


Table 4. New Specific Harmonics of the Inner Race Bearing Fault in Stator Current (Level 1/3 and 2/3) 


AN 1 
Sintece = (= (l-s)- if = kfin T P Í. 


Level 1/3 Level 3/3 

aA es 159.28 Hz 162.3 Hz 
le 259.28 Hz 262.3 Hz 
ayo 198.48 Hz 202.83 Hz 
ee 298.48 Hz 302.83Hz 


The motor is tested again under the same conditions monitoring for healthy case (load speed 1420 
Tr/min), with increasing the defect angle to level (3/3). 

Note that despite the decrease of the load from the level (1/3), the inner race bearing fault harmonics 
do not lose their amplitudes and remain visible, which shows the influence of the increase angle fault 
confirming that which was predicted in section III. Some examples show this difference in Figures 13 (a), 
(b), (c) and Figures 14 (a), (b), (c). 


Figure 13. (a)Experimental FFT spectrum of the stator current under Healthy condition (Black) s=0.053, 
faulty condition Level (1/3) (Blue)s=0.063, and faulty condition Level (3/3) (red)s=0.053 between 140Hz 
and 220Hz. (Zoom on (b) 159.28 Hz- 198.48 Hz Level (1/3) and (c) 162. 3 Hz-202.63 Hz Level (3/3)) 
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Figure 14. (a) Experimental FFT spectrum of the stator current under Healthy condition (Black), faulty 
condition Level (1/3) (Blue)s=0.063, and faulty condition Level (3/3) (red)s=0.053 between 250Hz and 
330Hz (Zoom on (b) 259.28 Hz- 298.48 Hz Level (1/3) and (c) 262. 3 Hz-302.63 Hz Level (3/3)) 


This two pairs of harmonics are summarized in Table 5. 


6. 


Table 5. New Specific Harmonics of the Inner Race Bearing Fault in Stator Current (Level 1/3 and 3/3) 


faee =((AN/ pX- 58) -1)F, — kf -( pf. 


Level 1/3 Level 2/3 
A=2, k=6 159.28 Hz 162.73 Hz 
259.28 Hz 262.73 Hz 


CONCLUSION 
A new comprehensive method for the calculation of different inductance values of an induction 


motor based on winding function approaches is presented. Inductances of healthy machine as well as 
dynamic eccentricity and inner race bearing fault is calculated and studied. The effect of the defect angle is 
verified. It is also shown that due to the air-gap asymmetry such as dynamic eccentricity or the inner raceway 
defect, the values of inductance coefficients have considerable variations with respect to the healthy 
condition and this last was applied for the study and detection of this fault. 


The experimental results which show current spectra of an induction machine with inner race 


defects, and are used to verify the new expression of the characteristic frequency of this fault. The 
experimental results clearly illustrate that the Stator Current Signature Analysis can be used to identify the 
presence of a bearing fault. 
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